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Abstract. We show that with a Next Generation Large Telescope one can 
detect the accelerated motions of ~ 100 stars orbiting the massive black hole at 
the Galactic center. The positions and velocities of these stars will be measured 
to astrometric and spectroscopic precision several times better than currently 
attainable enabling detailed measurements of the gravitational potential in the 
neighborhood of the massive black hole. We show that the monitoring of stellar 
motions with such a telescopes enables: (1) a measurement of the Galactic center 
distance i?o to better than 0.1% accuracy, (2) a measurement of the extended 
matter distribution near the black hole, including that of the exotic dark matter, 
(3) a detection of general relativistic effects due to the black hole including the 
prograde precession of stars and possibly the black hole spin, and (4) a detection 
of gravitational encounters between monitored stars and stellar remnants that 
accumulate near the Galactic center. Such encounters probe the mass function 
of the remnants. 


1. Introduction 

Near-infrared monitoring of stellar sources in the inner l" of the Galaxy with 
speckle imaging and adaptive optics techniques has enabled complete orbital 
reconstruction of se veral stellar somces orbiting the 4 x black hole at 

the Galactic center (ISchddel et al.l[2?i?i^. 1200,31: iGhez et al.ll2nf)Hh . Sources have 
been monitored with astrometric errors of a few milli-arcseconds, and radial 
velocity errors < 50 km s~^ (Eisenhauer et al. 2003; Ghez et al. 2003a), allowing 
the detection of the accelerated proper motions of ~ 10 stars. The star SO-2 has 
the shortest orbital period of ~ 15 yr. With current orbital data, the distance 
to the Galactic center Rn is mea sured to an accuracy of k, 5% (lEisenhauer et all 
120031 . 1200, 5l: see ISalim fc Gouldl[l99fll for a discussion of the physical principles 
underlying the measurement of Rq from the orbital data). 

Although the stellar orbits have provided unequivocal proof of the exis¬ 
tence of a massive black hole at the Galactic center, the matter content in the 
vicinity of the black hole remains largely unknown. Due to the short radial 
diffusion time as compared with the age of the bulge, a large number of massive 
compact remnants (5 —IOMq black holes ) could have segregated into the dynam- 
ical sphere of influence of the black hole (|Morri ^ Il99,3l: iMiralda-Escude fc Gouldl 
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2nnnt ). Furthermore, adiabatic growth of the massive black hole could have com¬ 


pressed a pre-existing distribution of cold dark matter (CDM) JE ser fc Sikivi'3 
Il987l: Quinlan, Hern auist , &: Sigurdssonifl995l: iGondolo fc Silldll999t) and stars 
(iPeeble^ll 97^ lYounelll98(ll into a density “spike”, although a variety of dynami¬ 
cal processes are capable of destroying such a spike (Ullio. Zhao, fc Kamionkowski) 
120011: iMerritt et alJ[2?)?i^: iGnedin fc Prima^l2004 : lMerrittll2004l . IfThe GDM 
consists of weakly-interacting massive particles (WIMPs), a sustained CDM 
spike might produce detectable WIMP annihilation radiation. 

We show here that the finer angular resolution and increased light-gathering 
power of a Next Generation Large Telescope (NGLT) will enable the detection 
of accelerated proper motions of ~ 100 stars with astrometric and spectro¬ 
scopic errors several times smaller than currently possible. By fitting to mock 
NGLT orbital data we determine the constraints such orbits place on the grav- 
itational potential near the massive bl ack hole at the Galactic center (see also 
Weinberg. Milosavlievic. &: Ghe3l2005l i. We demonstrate that the observations 
will yield a measurement of the density profile in a dark matter spike. They 
will also reveal the lowest order general relativistic effects and possibly higher 
order effects such as the black hole spin. We show that Rq will be measured to 
better than 0.1% precision, helping to place tight constraints on the shape of the 
Galactic dark matter halo. We also calculate the rate at which perturbations 
from massive remnants that accumulate near the Galactic center deflect the ob¬ 
served stellar orbital motions and describe how the remnants’ mass function can 
thereby be extracted from the monitoring data. We find that if stellar-mass 
black holes dominate the matter density then in 10 yr of monitoring at an as¬ 
trometric resolution of 0.5 mas approximately 10% of all monitored stars will 
experience detectable deflections in their orbital motions. 


2. Monitoring Stars near the Galactic Center with an NGLT 


2.1. Astrometric and Spectroscopic Limit 

The ability to extract information about the Galactic center environment via 
the monitoring of stellar motions is limited by the astrometric and spectroscopic 
precision of the measurements. We now estimate each in turn. 

Adaptive optics will enable an NGLT to operate near its diffraction limit. 
The astrometric limit achie vable with adaptive op tics on an NGLT scales with 
telescope aperture as (HE ao fc GolavitallT992^ . Since the astrometric limit 

of bright stellar sources at the Galactic center achieved today with 10 meter 
telescopes is 69iq ~ 1 mas, a 30 meter NGLT should achieve <5030 ~ 0.5 mas. 
This is a conservative estimate as future adaptive optic platforms may provide 
astrometric limits ~ 1% of the diffraction limit (M. Fitzgerald, private commu¬ 
nication), corresponding to 8630 ~ 0.1 mas in the A-band. 

Integral field spectroscopy with an NGLT is expected to provide a spectral 
resolution R = A/AA ~ 1 —2 x 10^. Thus, with an NGLT one can measure stellar 
radial velocities to an accuracy of Sv ~ 10 km s“^. This too is a conservative 
estimate. The faint, cool stars detectable with an NGLT will likely exhibit 
rich spectral features i ncluding mol e cular lines, enabling high spectral resolution 
studies. For example, iFiger et'aD (2003) measured the radial velocities of 85 
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cool stars in the central parsec of the Galaxy with velocity errors of ~ 1 km s ^ 
Thus, an NGLT may achieve 6 v ten times smaller than our above estimate. 


2.2. Confusion Limit 

The greatest obstacle to detecting and following hitherto unseen stars in the cen¬ 
tral 1" of the Galaxy is stellar crowding. The light of bright stars and underlying 
faint stars flood the pixel elements and impose a limit to the faintest detectable 
star. The limit depends on the tel escope specifications and on the lumino sity 
function of the stars in the field. In IWeinberg. Milosavlievic. fc Ghed (|2005l i we 
estimate the minimum luminosity a star must have in order to be detected and 
monitored with an NGLT by evaluating the confusion noise as a function of the 
stellar Lf- band luminosity func tion (KLF). We normalize the KLF to the obser¬ 
vations of ISchodel et al.l ((20031 1 of bright stars within the inner arcsecond and 
use the best-fit KLF slope from iGenzel et ahl (|20033). We integrate the KLF 
over stars brighter than the derived minimum luminosity to obtain an estimate 
of the total number of stars that can be detected and monitored with an NGLT. 
We find that a 30 meter NGLT will resolve N ~ 100 stars out to a magnitude 
limit of K 22 (Figure nj. For such an NGLT, confusion noise excludes the 
detection of stars with semi-major axes smaller than a ~ 300 AU and the ac¬ 
celeration threshold below which orbital motion is not detectable excludes stars 
with a > 3000 AU. 



Figure 1. The A-band magnitude limit and number N of stars with de¬ 
tectable orbital motions as a function of the aperture of a diffraction limited 
NGLT. Results are sho wn for power-law A-b and luminosity functions normal¬ 
ized to observa tions bvIScbp de] et alJ 1(200.' ) with slopes matching the ^ 2g 
range found bv IGenzel et al.l ll2003all . After Weinberg. Milosavlievic. fc Ghed 

(l20^ . 
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3. Modeling the Orbital Dynamics 

Our goal is to estimate the uncertainties in the parameters that describe the 
gravitational potential through which the monitored stars move. We do so by 
first generating mock orbital data assuming a particular model for the potential 
(e.g., black hole mass, extended matter density profile, etc.) and a particular 
observing program (e.g., time baseline of 10 yr, 10 observations per year per 
star with measurement errors [56, 6v] = [0.5 mas, 10 km s~^]). The orbital pa¬ 
rameters are drawn from the distribution function of a realistic power-law 

density profile assuming randomly oriented orbits and considering only those 
orbits detectable with an NGLT (e.g., N = 100 stars with semi-major axes in 
the range 300 AU < a < 3000 AU). 

3.1. Equations of Motion 

An orbit fitting routine must model various non-Keplerian motions in order to 
extract information about the gravitational potential near the massive black 
hole. Effects that give rise to deviations from Keplerian motions include the 
Newtonian retrograde precession of an orbit due to the presence of an extended 
matter distribution, the relativistic prograde precession, and the frame drag¬ 
ging effects due to the black hole spin. The post-Newtonian approximation to 
the geodesic equation for test particles orbiting a spherically symmetric mass 
accounts for the above effects and provides an appropriate d escription of the 
observable motions fsee lWeinberg. MilosavlieviC' fc Ghezlliin^ . 

We also examine the prospect of directly measuring the mass function of 
the remnants by looking for changes in orbital semi-major axis due to the gravi¬ 
tational interaction between monitored stars and background massive remnants. 
Such interactions will be of two types: strong encounters and weak encounters. 
The former correspond to cases where the impact parameter of an encounter sat¬ 
isfies d < where m* is the mass of the background star and Urei is the 

relative velocity of the encounter, causing a large change in semi-majo r axis. We 
calcul ate the rate of strong encounters by extending the analysis of iGoodmanI 
( 19831 ). who considered the encounter rate in globular clusters. Weak encoun¬ 
ters describe the slow diffusion of a star’s orbital energy due to the continual 
small angle scattering by the background stellar distribution, whose total po¬ 
tential is the sum of disc rete point-mass potentia ls. We solve the orbit-averaged 
Fokker-Planck equation (|Gohn fc Kulsru^ll978l i to calculate the rate at which 
weak encounters change a star’s semi-major axis. The rate of strong and weak 
encounters is proportional to the mass of the remnants. 

3.2. Parameter Estimation 

A star’s projected orbit is described by six phase space parameters. The poten¬ 
tial contributes six additional parameters: the black hole mass, its 3-dimensional 
position, and the normalization and slope of the extended matter distribution. 
The dimension of the parameter space of our model is therefore J = QN -\- 6. 
Thus, to extract parameter constraints from the N ~ 100 stars detectable with 
a 30 meter NGLT requires exploring a J ~ 600 dimensional parameter space. 
Since the computational cost of the grid-based approach increases exponentially 
with J, parameter estimation on a J-dimensional grid is not practical for N ^ 3. 
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Instead we use the Markov Chain Monte Carlo (MCMC) method which provides 
an efficient means of mapping out the probability distribution of large dimen- 
si onal models. The details of our implemen tation of the MCMC method is given 
in IWeinberg. Milosavlievic. fc GhezI (l2nn,^ . 


4. Results 

We show how well an NGLT constrains the mass and distance to the massive 
black at the Galactic center and the extended matter distribution of stars and 
dark matter around the black hole. We also demonstrate that relativistic effects 
are detectable with an NGLT. Finally, we show that interstellar interactions 
between monitored stars and background massive remnants will be observed 
with an NGLT if the population of stellar-mass black holes predicted by theory 
exists. Although we show results for a diffraction limited 30 meter NGLT with 
{S6,6v) = (0.5 mas, 10 km s~^) and {66, 5v) = (0.1 mas, 1 km s~^), since the 
parameter uncertainties scale with measurement error a (i.e., 69 and 6 v) and 
number of monitored stars N as the results can be used to describe 

the capabilities of an NGLT with different specifications. For example, based 
on Figure ^ a 100 meter NGLT will detect ten times as many stars as a 30 
meter NGLT. Thus, if a 100 meter telescope has astrometric and spectroscopic 
errors that are smaller than those of a 30 meter telescope by a factor of five, the 
parameter uncertainties will be ~ lOx smaller. 

4.1. Measuring Mbh and Rq 

For an astrometric limit of 69 = 0.5 mas and a spectroscopic limit of 6 v = 
10 km s~^ the fractional uncertainties in Mbh and Rq are less than 0.1% at the 
99.7% level (FigurejJl), a factor of ~ 100 times better than present uncertainties. 
For astrometric and spectroscopic limits that are a factor of five smaller the 
fractional uncerta inties in Mbh and R n are smaller by almost a factor of five. 

As shown by lOlling fc Merrifieldl ()2000l . |20^) measuring the Galactic con¬ 
stants Rq and the Galactic rotation speed ©o to high accuracy will constrain 
the shortest-to-longest axis ratio q = c/a of the Galactic dark matter halo to 
similar accuracy. The shape parameter q is an important diagnostic of dark 
matter models and structure formation and is currently poorly constrained in 
all galaxies including the Milky Way. 

4.2. Measuring the Extended Matter Distribution 

In Figure 121 we show the constraints an NGLT will place on the extended matter 
distribution for input power-law models with M^^i{r < 0.01 pc) = OOOOMq 
and 7 = 1.5 or 7 = 2. We find that an NGLT will measure Mext and 7 to 
20 — 30% accuracy. We find that an extended matter distribution is detectable 
(i.e., observations yield a lower bound) for 69 = 0.5 mas and 6 v = 10 km s~^ 
as long as Mext(?’ < 0.01 pc) ^ 15OOM0. Thus, if the dark matter distribution 
matches theoretical expectations and forms a density spike, its influence on 
the orbits will be detectable with an NGLT. A detection would constitute a 
measurement of the gravitational influence of dark matter on the smallest scales 
yet. 
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Figure 2. Left panel: The constraints on Mbh and Rq obtainable with an 
NGLT. The solid contours show the 68%, 95%, and 99.7% confidence levels 
assuming an astrometric limit of 60 = 0.5 mas and a spectroscopic limit of 
Sv = 10 km s~^ for a mock sample of stars detectable with an NGLT. The 
line-dot contour shows the 99.7% confidence level for a different draw of mock 
orbits. The dashed contour shows the 99.7% confidence level for smaller 
astrometric and spectroscopic limits of 60 = 0.1 mas and 6v = 2 km s~^. 
Right panel: The constraint on the extended matter distribution obtainable 
with an NGLT. The input models have power-law slope of 7 = 1.5 and 7 = 2 
and an input enclosed mass of 600 0Mrb within 0.01 pc. The line styles ar e 
the same as in the left panel. After [Weinberg. Milosavlievic. fc “Gh^ll^OOl . 


4.3. Measuring Relativistic Effects 

The relativistic prograde precession (RPP) causes a pericenter advance per rev¬ 
olution of A(/>pro = SirRx/a d — e^), where R .g = 2GMbh/c^ is the Schwarzschild 
radius of the black hole (see IWeinberdHoT^ . The apparent apocenter shift per 
revolution caused by the RPP is As A(/)proa(l -|- e)/Ro = SttRs/R od — e), 
which corresponds to an apparent shift of ~ 1 mas for the star SO-2. Consider 
an orbit seen face on and observed for Aorb complete periods. Since the pre¬ 
cession angles per revolution add linearly, the signal-to-noise from the RPP is 
Spro ~ AsNorh/^0, or 


Spro ~ 


0.1 


Aorb / Mbh \ / Rq 
1 — e V4x 10® Mq / V 8 kpc 


-1 


66 


1 mas 


-1 


( 1 ) 


If only one star with e > 0.96 is monitored over a single period at an astrometric 
accuracy 69 = 0.5 mas, the RPP shift will be measured to 5 -cj accuracy. Several 
such stars are expected to be detected in a 30 meter NGLT’s sample of 100 stars. 

Frame dragging effects due to a spinning black hole also cause a precession 
of a stellar orbit’s pericenter. The spin precession per revolution for a star 
orbiting a black hole with spin angular momentum J is given approximately by 
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(see IWeinberdHoT^ equation (9.5.22); note different notation) 


^'/’spin 


-Svrj 


/ GMbh 
V cL 


3 

COS Ip 


( 2 ) 


where ip is the angle between the orbital angular momentum vector and the 
black hole spin axis and 0 < j = cJ/GM^^ < 1 is the black hole spin parameter. 
The black hole spin induces an apocenter shift that is smaller than the RPP 
shift by a factor of ~ u/c. For an orbit observed face-on the signal-to-noise from 
a spin-induced apocenter shift is 


-'Spin 


^ jWorb cosip f Mbh 

V'(l + e)(l -e)a U X 106 Mq) 

X ( ^ / Rq / se 

VlOOOAU/ VSkpc/ Vlnias/ 


( 3 ) 


In order to be detectable a spin-induced orbital precession requires an astromet¬ 
ric precession 66 ^ 0.05 mas or th e favorable detection of a star on a highly co m- 
pact and eccentric orbit (see aLso l.Taroszvnskil[l998l:l^aEfile fc Mathewsl|200flr) . 

To confirm the above order of magnitude estimates we generated mock 
orbits using the post-Newtonian corrections to the equations of motion and fit 
to these orbits using the MCMC method while treating the speed of light as 
a parameter in the model. We found that the RPP effect is measurable for 
69 = 0.5 mas (speed of light constrained to 5% accuracy) while the spin induced 
precession is not. 



Figure 3. The probability that after T = 10 yr the semi-major axis of 
a star is changed from an initial value ai to a final value a/ due to strong 
and weak encounters. The three panels show results for the cases ai = 200, 
1000, 2500 AU with e = 0 and e = 0.9. The dotted lines show the separate 
contributions from strong and weak encounters and the solid lines show the 
net probability. The vertical dashed lines are the uncertainties in a, assuming 
66 = 0.5 mas. 


4.4. Interstellar Interactions 

Milosavlicvic fc Wcinbei^ (|2006l i estimate the rate at which monitored stars 
experience detectable encounters with background stars and stellar remnants 
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(see also I Weinberg. Milosavlievic. fc Ghedl2nn5l l. Such encounters will manifest 
themselves as changes in a star’s orbital energy, i.e., semi-major axis (Figure 
0, at a rate proportional to the mass of the background sources. We found 
that if the background sources are dominated by stellar-mass black holes, as 
predicted by estimates of mass segregation in the vicinity of a massive black hole 
( Morrisl[T99,‘ll: iMiralda-Escude fc Gouldil20nnll , approximately 30 encounters will 
be detected over a ten year baseline for 66 = 0.5 mas. 


5. Conclusions 

We have shown that the monitoring of stellar orbits around the massive black 
hole at the Galactic center at the high astrometric and spectroscopic resolution 
attainable with an NGLT enables one to probe the deep gravitational potential of 
the region. Many exciting measurements are achievable even for modest (factor 
of a few) improvements over the astrometric capabilities of current 10 meter 
class telescopes. The future success of Galactic center research greatly depends 
on advances in the astrometric capabilities of the next generation of telescopes. 
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